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Abstract 

We emphasize the sizeable effects of absorption on high-energy 'soft' processes, and, 
hence, the necessity to include multi-Pomeron-Pomeron interactions in the usual multi- 
channel eikonal description. We present a model which includes a complete set of the 
multi-Pomeron vertices and which accounts for the diffusion in both, the impact parameter 
and ln(kt), of the parton during its evolution in rapidity. We tune the model to the 
available data for soft processes in the CERN-ISR to Tevatron energy range. We make 
predictions for 'soft' observables at the LHC. 



1 Motivation 

There are three main reasons for revisiting soft pp high energy interactions at this time. 

A. This paper is concerned with the description of the high energy behaviour of "soft" 
observables such as Utot, dae\/dt, dasB / dtdM"^ , particle multiplicities etc. in terms of basic 
physics. This physics predated QCD and is sometimes regarded as the Dark Age of strong 
interactions. However, it is unfair to call this the Dark Age. We had a successful description 
of these processes in terms of the exchange of Regge trajectories linked to particle states in 



the crossed channels [1] . The dominant exchange at high energy is the Pomeron, and we have 
Gribov's Reggeon calculus [2] to account for the multi-Pomeron contributions. However the 
available data did not reach high enough energy to distinguish between the different scenarios 
[3, 4] for the high-energy behaviour of the interaction amplitude [5, 6]. 

In the 'weak coupling' scenario the total cross section atot{s — > oo) — > const, and in order 
not to violate unitarity, and to satisfy the inequality 



the triple- Pomeron vertex must vanish when t — > 0, that is the triple- Pomeron coupling g^p oc t. 
In this case, the large logarithm coming from the integration over the mass of the system 
produced in diffractive dissociation (/ dM^/M^ ~ Ins) is compensated by the small value of 
the mean momentum transferred through the Pomeron, (t) oc 1/lns. 

On the other hand, in the 'strong coupling' scenario, where Utot (Ins)'' with < 77 < 2, 
the inequality (1) is provided by a small value of the rapidity gap survival factor which 
decreases with energy. 

The present diffractive data arc better described within the 'strong coupling' approach [7], 
and in this paper we shall give predictions for the LHC for this scenario. However, the possibility 
of the 'weak coupling' scenario is not completely excluded yet. Therefore, it is quite important 
to study the different channels of diffractive dissociation at the LHC in order to reach a final 
conclusion and to fix the parameters of the model for high-energy soft interactions. So the 
first motivation is the intrinsic interest in obtaining a reliable, self-consistent model for soft 
interactions, which may be illuminated by data from the LHC [8] . 

B. In turn, obtaining a reliable model will be of great value for predictions of the gross 
features of soft interactions. In particular, it is essential for understanding the structure of the 
underlying events at the LHC. 

C. The third reason for studying soft interactions arises because it may not be an easy task 
to identify the production of a new object at the LHC when it is accompanied by hundreds of 

other particles emitted in the same event. For the detailed study of the new object. A, it may 
be better to select the few, very clean, events with the Large Rapidity Gaps (LRG) on cither 
side of the new object, see, for example, [9, 10, 11, 12]. That is to observe the exclusive process 
pp — > p + A+p. In such a Central Exclusive Process (CEP) the mass of A can be measured with 
very good accuracy {AMa ~ 1 — 2 GeV) by the missing-mass method by detecting the outgoing 
very forward protons. Moreover, a specific — selection rule [13] reduces the background 
and also greatly simphfies the spin-parity analysis of A. However, the CEP cross section is 
strongly suppressed by the small survival factor, 5"^ <^ 1, of the rapidity gaps. Thus we need a 
reliable model of soft interactions to evaluate the corresponding value of [14]. Moreover, it 
is important to have a model which contains t-channel components of different size in order to 
evaluate the possible effects of the 'soft-hard factorisation' breaking. This is the subject of the 
following paper [15]. 
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Figure 1: The eikonal model of elastic scattering 



2 Resume of the eikonal formalism 
2.1 Single-channel eikonal model 

First, we briefly recall the relevant features of the single- channel eikonal model. That is, 
we focus on elastic unitarity. At high energy the position of the fast particle in the impact 
parameter, b, plane is to a good approximation frozen during the interaction, since the value 
of b is fixed by the orbital angular momentum I = 6-\/i/2 of the incoming hadron. There is no 
mixture between the partial wave amplitudes with different /. The well known solution of the 
elastic unitarity equation, 

2ImTei(s,6) = \T,,{s,b)\' + G,Us,b), (2) 

may be written in terms of the phase shift Si as 

Sci = l + iTei = e^'^', that is T^i = i{l - e^'^') . (3) 

The presence of inelastic channels, given by Ginei in (2), leads to the phase Si having an imag- 
inary part. That is. Si becomes a complex number. Moreover, at high energies we know that 
ReTei/ImTei is small. 

Now, in the framework of the eikonal model, the elastic amplitude, 

Tel = ^(1 - (4) 

is obtained by the sum of Regge-exchange diagrams, which is equivalent to the iteration of the 
elastic unitarity equation, (2), as shown in Fig. 1. In other words, s-channel elastic unitarity 
gives 

lm%,{s,b) = l-e-^/2 (5) 
= (l-e-^/2)^ (6) 
a,Us,b) = 1-e-^, (7) 

where Q{s, 6) > is called the opacity (optical density) or eikonal^; Q/2 plays the role of —2iSi. 
It is the Fourier transform of the two-particle (s-channel) irreducible amplitude, A{s,qt). That 
is^ 

^is,b) = ^ j d^qtA{s,qt)e^^-'" , (8) 



^Sometimes £7/2 is called the eikonal. 

^We use the bold face symbols qt and b to denote vectors in the transverse plane. 
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where — —t, and where the amphtude is normahzed by the relation a"tot(s) = ImTei(s, t — 0). 
From (7), we sec that exp(— is the probabihty that no inelastic scattering occurs at 
impact parameter b. 

After summing over b (that is, all partial waves) we obtain the total, elastic and inelastic 
cross sections 

(Jtot = 2 J d''blmn,is,b) ^ 2jd% (l-e-"/2) (9) 
aei = I d%\T,,{s,b)\' = I d% {1 - e-^'/'Y (10) 
(Tinei = J d%[2ImT4s,b)-\T4s,b)\^] = jd%{l-e-''). (11) 

Below we neglect the imaginary part of fi, apart from the contribution of secondary Reggeons 
to high-mass diffractive dissociation. At high energies, the ratio ReTei/ImTei is small, and can 
be evaluated via a dispersion relation. 



2.2 Inclusion of low- mass diffractive dissociation 

So much for elastic diffraction. Now we turn to inelastic diffraction, which is a consequence 

of the internal structure of hadrons. Besides the pm^c elastic two-particle intermediate states 
shown in Fig. 1, there is the possibility of proton excitation, p — >■ N*. As a rule such excitations 
are not included in the opacity ^2, but are treated separately 

This is simplest to describe at high energies, where the lifetime of the fluctuations of the 
fast proton is large, r ~ Ejrr?^ and during these time intervals the corresponding Fock states 
can be considered as 'frozen'. Each constituent of the proton can undergo scattering and thus 
destroy the coherence of the fluctuations. As a consequence, the outgoing superposition of 
states will be different from the incident particle, and will most likely contain multiparticle 
states, so we will have inelastic, as well as elastic, diffraction. 

To discuss inelastic diffraction, it is convenient to follow Good and Walker [16], and to 
introduce states (pk which diagonalize the T matrix. Such eigenstates only undergo elastic 
scattering. Since there are no off-diagonal transitions, 

(0,|T|0fc)=O iorj^k, (12) 

a state k cannot diffractively dissociate into a state j. We have noted that this is not, in general, 
true for hadronic states, which are not eigenstates of the ^-matrix, that is of T. To account for 
the internal structure of the hadronic states, we have to enlarge the set of intermediate states, 
from just the single elastic channel, and to introduce a multichannel eikonal. We will consider 
such an example below, but first let us express the cross section in terms of the probability 
amplitudes of the hadronic process proceeding via the various diffractive eigenstates^ (pk- 

^The exponent exp(— fife) describes the probability that the diffractive eigenstate is not absorbed in the 
interaction. Later we will see that the rapidity gap survival factors, S"^, can be described in terms of such 
eikonal exponents. 



3 



Let us denote the orthogonal matrix which diagonahzes ImT by a, so that 



ImT = aT^J with (0j|T^|0fc) = Sju- (13) 
Now consider the diffractive dissociation of an arbitrary incoming state 

fc 

The elastic scattering amplitude for this state satisfies 

O'lImTli) = Y.M'T^ = (T^^^), (15) 

where = {(t)k\T^\(t>k) and where the brackets of {T^) mean that we take the average of 
over the initial probability distribution of diffractive eigenstates. After the diffractive scattering 
described by Tfj, the final state |/) will, in general, be a different superposition of eigenstates 
from that of which was shown in (14). At high energies we may neglect the real parts of 
the diffractive amplitudes. Then, for cross sections at a given impact parameter 6, we have 

- 2ImO-|Tb) = 2 E M'T^ = 2(T^) 



d?h 

C^CTel + SD 



It follows that the cross section for the single diffractive dissociation of a proton. 



d?h 

is given by the statistical dispersion in the absorption probabilities of the diffractive eigenstates. 
Here the average is taken over the components k of the incoming proton which dissociates. If 
the averages are taken over the components of both of the incoming particles, then in (17) we 
must introduce a second index on T^, that is T^, and sum over k and i. In this case the sum 
is the cross section for single and double dissociation. 

At first sight, enlarging the number of eigenstates we may include even high- mass proton 
dissociation. However here we face the problem of double counting when the partons originating 
from dissociation of the beam and 'target' initial protons overlap in rapidities. For this reason 
high-mass dissociation is usually described by "enhanced" multi-Pomeron diagrams. The first 
and simplest is the triple-Pomeron graph, see Fig. 5 below. 
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3 Triple-Regge analysis accounting for absorptive effects 



The total and elastic proton-proton cross sections are usually described in terms of an eikonal 
model, which automatically satisfies s-channel elastic unitarity. To account for the possibility 
of excitation of the initial proton, that is for two-particle intermediate states with the proton 
replaced by N*, we use the Good- Walker formalism [16]. Already at Tcvatron energies the ab- 
sorptive correction to the elastic amplitude, due to elastic eikonal rescattering, is not negligible; 
it is about —20% in comparison with the simple one Pomeron exchange. After accounting for 
low-mass proton excitations (that is A^*'s in the intermediate states) the correction becomes 
twice larger (that is, about —40%). Indeed, the possibihty of proton excitation means that we 
have to include additional inelastic channels which were not accounted for in the irredicible 
amplitude A of (8). This enlarges the probability of absorption for the elastic channel, that is 
the effective opacity fl. In terms of the Good- Walker formalism, the stronger absorption follows 
from the inequality^ 

(Ofee-^*) < (Ofc)(e-^'=), (18) 

where we average over the diffractive (Good- Walker) eigenstates. 

Next, in order to describe high-mass diffractive dissociation, daso/dM'^, we have to include 
an extra factor of 2 from the AGK cutting rules [17]. Thus, the absorptive effects in the triple- 
Regge domain arc expected to be quite large. The previous triple-Regge analyses (see, for 
example, [18]) did not allow for absorptive corrections and the resulting triple-Regge couplings 
must be regarded, not as bare vertices, but as effective coiiplings embodying the absorptive 
effects [19]. Since the inelastic cross section (and, therefore, the absorptive corrections) expected 
at the LHC are more than twice as large as that observed at fixed-target and CERN-ISR 
energies, the old triple-Regge vertices cannot be used to predict the diffractive cross sections 
at the LHC. 

Thus, it is necessary to perform a new triple-Regge analysis that includes the absorptive 
effects explicitly. Such an analysis has recently been performed [7] using the fixed-target FNAL, 
CERN-ISR and Tevatron data that are available in the triple-Regge region. The 'PPP', 'PPR\ 
^RRP\ ^RRK' and ttttP contributions, were included assuming either the 'strong' or 'weak' 
coupling scenarios for the behaviour of the triple- Pomeron vertex. To account for the absorptive 
corrections a two-channel (Good- Walker) eikonal model was used, which describes well the total, 
(Jtot, and elastic, da^i/dt, pp and pp cross sections. 

In the 'strong' coupling case, a good x^/DoF=167/(210-8)=0.83 was obtained. In compar- 
ison with the old triple-Regge analysis [18], a twice larger relative contribution of the 'PPi?' 
term was found. This is mainly due to the inclusion of the higher-energy Tevatron data in the 
analysis. 

''When we go from a single- to a many-channel eikonal, we may write = (17) +5^ with {5^) ~ 0. It follows 
that (fifce~^'=) = (r2)(e~^'=) -|- (Jfee~*'=)e~^^^ which, since the second term is negative, is less than (ri)(e~^'=). 
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Since the absorptive effects are included explicitly, the extracted values of the triple- Reggeon 
vertices are now much closer to the bare triple-Regge couplings. In particular, the value 



gppp = Xqn, 



where A ~ 0.2 



(19) 



is consistent with a reasonable extrapolation of the perturbative BFKL Pomeron vertex to the 
low scale region [20]; here qn is the Pomeron-proton coupling. Note also that these values of 
the ^PPP^ and ^PPR vertices allow a good description of the HERA data [21] on inelastic 
J/V' photoproduction, ^ J/ip + Y , where the screening corrections are rather small. 

The 'weak' coupling scenario leads to a larger x^/DoF=1.4 and to a worse description of 
the 7P ^ J/ip -\-Y process at the lowest values of t. At the LHC energy the 'weak' coupling 
fit predicts about 3 times smaller inclusive cross section da^^/dtdM^ at ^ = M'^/s = 0.01 and 
low t in comparison with that predicted in the 'strong' coupling case. 

4 Model with a complete set of multi-Pomeron vertices 

Note that the effects due to the triple- Pomeron vertex (19) are rather large. Indeed, the contri- 
bution caused by such vertices is enhanced by the logarithmically large phase space available in 
rapidity. In particular, the total cross section of high-mass dissociation is roughly^ of the form 



where A reflects the suppression of high-mass dissociation in comparison with elastic scattering 
and the Ins factor comes from the integration / dM^/M^ ~ Ins. Thus actually we deal with 
the parameter Alns > 1 at collider energies. For each fixed rapidity interval the probability of 
high-mass dissociation (or, in other words, the contribution due to the triple-Pomeron vertex) 
is relatively small. However the cumulative effect in the complete interaction amplitude is 
enhanced by the large phase space available in rapidity. 

As a consequence, the contribution of the corresponding, so-called 'enhanced', diagrams, 
with a few vertices, is not negligible. Moreover, we cannot expect that more complicated 
multi-Pomeron interactions, driven by the gll^ vertices, which describe the transition of n to m 
Pomerons of Fig. 2, will not affect the final result. It looks more reasonable to assume that 
gf^ oc X^+™- than to assume that g'^ = for any n + m > 3. Thus we need a model which 
accounts for the possibility of multi-Pomeron interactions (with arbitrary n and m). 

In this paper we extend and develop the partonic approach of Ref . [22] . While the eikonal 
formalism describes the rescattering of the incoming fast particles, the enhanced multi-Pomeron 
diagrams represent the rescattering of the intermediate partons in the ladder (Feynman dia- 
gram) which describes the Pomeron-exchange amplitude. 

^Here, for simplicity, we assume an essentially flat energy dependence, a ~ with elns < 1. 




(20) 
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Figure 2: A multi-Pomeron vertex 





n{y = 0) 

Figure 3: The evolution of the elastic bare Pomeron amplitude, Qkiy,b). 



Indeed, we start with the generation of ladder-type structure of 'elastic' bare Pomeron 
exchange amplitude. It may be generated by the evolution equation [23] (in rapidity, y, space) 



dn{y,b) 
dy 



(21) 



where b is the two-dimensional vector in impact parameter space. A is the probability to 
emit new intermediate partons (denoted c) within unit rapidity interval; it is analogous to the 
splitting function of DGLAP evolution. The impact parameter of parton c is not frozen in the 
evolution. At each step b can be changed by a constant amount Ab in any direction, leading 
to the diffusion represented by the second term on the right-hand side of (21) where a' plays 
the role of the diffusion coefficient [24]. The evolution is shown symbolically in Fig. 3. The 
solution of (21) is 

n{y, b) = no exp{yA - b"" / Aa'y) / Ana'y. (22) 

It represents the opacity (at point y, 6), corresponding to the incoming particle placed at 6 = 
and y = 0. 

It may be helpful to explain why (21) was written in terms of the opacity Vt. First we 
note that the discontinuity of the amplitude generated by (21) does not contain a two-particle 
s-channel intermediate state; it corresponds to a pure inelastic high multiplicity process, see 
Fig. 4. Due to elastic unitarity, (2) this inelastic interaction leads to elastic pp scattering. Thus 
we have to put the solution of (21) into the eikonal formulae of (5)-(7), as it does not contain 
two-particle s-channel states. 
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Figure 4: A pure inelastic high multiphcity process 
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Figure 5: The ladder structure of the triple-Pomeron amplitude. 

In momentum space the solution (22) corresponds to the amplitude 

A{s,t) = Ao 5^+^+^^'*. 



(23) 



That is to the bare Pomeron exchange amplitude, where the Pomeron trajectory has intercept 
a(0) = 1 + A and slope a'. 

A multi- Pomeron enhanced contribution arises from the absorption of intermediate s-channel 
partons c during the evolution of Q in y. The simplest example is the triple-Pomeron diagram 
in which parton c undergoes an extra rescattering with the target parton k, as shown in Fig. 5. 
Allowing for many rescatterings, we have to sum over different numbers of ladders between 
partons c and k. Assuming an eikonal form for the multi-Pomeron-proton vertex, it is natural 
to replace (21) by 



dy 



(24) 



where the 'opacity' Qk describes the transparency of the target k. As we are dealing with the 
elastic amplitude we use e~^^*/^ and not e'^^'' . The coefficient A arises since parton c will have 
a different absorption cross section from that of eigenstate i. Naively, we may assume that the 
beam i contains a number 1/A of partons. The factor e~'^^'=/^ generates multi-Pomeron vertices 
of the form 

= n-m- A"+"^2^iv/2 for n + m > 3 . (25) 

where g]^ is defined in Fig. 2. Even though A ^ 0.25, the role of the factor e"'^^''^'^ is not 
neghgible, since the suppression effect is accumulated throughout the evolution. For instance. 



8 



if A <S 1 the full absorptive correction is given by the product AQy/2, where the small value of 
A is compensated by the large rapidity interval Y. 

In terms of Regge diagrams, (24) sums up the system of fan diagrams in which any number 
m of "lower" Pomerons couples to a fan vertex g^^. Such multi-Pomeron diagrams are called 
"enhanced" , since their contribution is enhanced, in comparison with the eikonal diagrams, by 
the large available rapidity interval Y. In order to include the rescattering with the beam i we 
replace (24) by 



The final term in the exponent is the opacity of the beam i, which depends on the rapidity 
interval y' — Y — y, with Y — Ins. The equation for the opacity Qi has an analogous form 



in which we now evolve in the opposite direction starting from the boundary condition Q(y' = 0) 
at y = y. 

Recall that the fit to the data in the triple- Reggeon domain indicated a very small (consistent 
with zero) t-slope of all the triple- Reggeon vertices [7, 18]. Thus, as the size of the multi-Rcggcon 
vertices are negligible in comparison with the size of the incoming hadron, we may write the 
absorptive corrections (that is, the exponential factors on the right-hand-side of (26,27)) such 
that the opacities are taken at the same point in the impact parameter plane b. 

Since the intermediate parton may be absorbed by the interaction with the particles (par- 
tons) from the wave function of both the beam or target hadron, we now need to solve the two 
equations, (26) and (27). This is done itcrativcly. Moreover, note that the opacities Q now 
depend on two vectors in impact parameter space - the separation bi between the position of 
the intermediate parton c and the beam hadron, and the separation b2 between c and the target 
hadron. The argument b in (26,27) now symbolically denotes both bi and b2. The resulting 
solution Q(y, b) is then used in the eikonal formulae for the elastic amplitude, giving 



A more detailed description of the amphtude, and the cross sections of the different diffractive 
processes can be found in Ref. [22]. 

4.1 Multi-components in both the s- and t-channels 

As mentioned above, as in [22], we use three diffractive components in the s-channel. In other 
words, we use a 3-channel eikonal for the rescattering of fast particles. The transverse size 
squared of each eigenstate is proportional to the corresponding absorptive cross section; i?f oc 
CTj. That is, we assume that the parton density at the origin is the same for each eigenstate. The 




(26) 




(27) 



Tei(b) = 1 - exp(-f](b = bi - b2)/2). 



(28) 
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shape of the Pomeron-nucleon vertex is parametrised by the form factor V{t) — e'^'^*' /{1 — t/ diY , 
whose Fourier transform, T^(b), plays the role of the initial conditions for Q(y = 0, b). However, 
now, we allow for a non-zero slope {a' ^ 0) of the (bare) Pomeron trajectory. 

A major development, of the model of [22], is that we use four different i-channel states, 
which we label a: one for the secondary Reggeon (i?) trajectory and three Pomeron states 
[Pi, P2, P3) to mimic the BFKL diffusion in the logarithm of parton transverse momentum, 
ln(A;t) [25]. To be precise, since the BFKL Pomeron [26] is not a pole in the complex j-plane, 
but a branch cut, we approximate the cut by three t-channel states of a different size. The 
typical values of kt in each of the three states is about kn ~ 0.5 GeV, kt2 ~ 1-5 GeV and 
kt3 ~ 5 GeV. Thus the system of evolution equations (26,27) is replaced by^ 



dntivMM) ^ ^-xinliyMM)+^ny'MM)]/2 ( ^ ^a^y^ b2) + Ka'^^^', (29) 

dy \ d^bi J 

"^^'^^ly^"^'^ ^ e-^[^^(^'^-^^)+^°(^'-^^'^^)l/2 (a'^ + a'^^^ ^I1{y\ b^, b^) + K^^Of , (30) 

where A" = q;(0) — 1 and a'^^ — a'p for a = Pi, P2, -P3, while for the secondary Reggeon, [a — R), 
which is built of quarks, we take A^ = ocr{Qi) — 0.6 and a'^ = 0.9 GeV~^, so that the last 
term VrrQ.^ is diagonal with Vrr = — 1 to account for the spin \ nature of quarks. The key 
parameters which drive the evolution are the intercepts A and the slopes a' . In general, each 
component a may have different values of A^ and a'^^. We discuss the values in Section 5. 

In the exponents, the opacities (Qfc) are actually the sum of the opacities (O^') with 
corresponding coefficients. Namely 

pp 

q'' = n^^v,, + n^v,,. (31) 



We chose Vpp = (1/3) since, at the leading order, the probability of interaction of two com- 

\2 



ponents of different size (kt) is proportional to the ratio (^42/^*1) - We take v'^^ = 1/27 since 
the third (smallest size) component collects all the higher kf contributions, and therefore here 
the mean value of kt is larger. For the screening of the secondary Reggeon by the Pomeron we 
take just the colour factor v^p = Cf/Ca = (4/9), as we assume that the secondary Reggeon 
is composed of a i-channel quark-antiquark pair. Finally the factors vrr — 1.8 and vrr — 4 



^Strictly speaking both opacities 0.% and Q.k depend on both subscripts i and k. Here we keep only one 
subscript to distinguish the parent hadron for each active parton (gluon). 
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were tuned to give a reasonable reproduction of the secondary Reggeon contributions to the 
available pp — > p + X data^. 

The transition factors Vaa' between the different i-channel components are fixed by the 
properties of the BFKL equation. The only non-zero factors, apart from Vrr = — 1, are 

KpiP2 = P^^Vpp, Vp^p^ = p^^Vpp, Vp^p, = p^%p, Vp^p^ = p^'^v'^p, (32) 

where 

p" is the density of partons emitted in the rapidity evolution of the i-channel component a. 
The remaining transition factors were set to zero. That is 

Vp,p, ^ Vpip, ^ 0, Viu^^VaR^O, and Ka' = for a = a' . (34) 

For each Good- Walker s-channel component (i.k). the initial conditions are fixed by the 
parton (matter) distribution in the corresponding diffractive eigenstate 

nt{y = 0, b) = = A y e^-^-^md\ , (35) 

where Pi{t) — 7j/3(7jt). The parametrisation P{t) — /^oe'^^V (l—t/di)'^ was used for the Pomeron, 
while for the secondary Reggeon we chose the Gaussian form (3{t) — Prc^'^^. 

The relative couplings (and the corresponding size) of the components were taken to (a) 
reproduce the cross section of low-mass dissociation measured at the GERN-ISR [27], and (b) to 
make all three components quite different from each other; 71 = 1.80, 72 = 0.82, 73 = 0.38. All 
the coefficients in the decomposition |p) = J2 Oil^i) are taken to be = l/-\/3 (with i = 1, 2, 3). 

To avoid possible double counting, and to fix the boundary between the 'low' and 'high' 
mass dissociation, we introduce a threshold Ay = 1.5 in rapidity for the actual start of the 
evolution of (29,30). That is we start the evolution at y = Ay = 1.5 and not aX y = 0. Hence 
the available rapidity interval becomes 6Y = ln(s) — 2Ay. Proton excitation (dissociation) 
which covers a rapidity interval larger than Ay (i.e. ln(M^/so) > 1-5) is called 'high-mass 
dissociation'. 

It is natural to separate the different contributions in terms of rapidity, since in QGD 
the interference between the different diagrams for gluon radiation leads to angular (rapidity) 
ordering of emitted gluons, at least to leading log accuracy^. 

Note that the initial condition (35) is only valid for the secondary Reggeon and for the large 
size Pomeron component {a — Pi). For the smaller size Pomeron components we use 

(y = 0, b) = (y = 0, b)^,, and Qf^(y = 0, b) = (y = 0, b)^,,<^. (36) 
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Figure 6: The irreducible amplitude Fik{Y,h) of a high energy interaction. 



4.2 Total and differential cross section formulae 

To calculate the elastic amplitude we need the s-channel two-particle irreducible amplitudes for 
the scattering of the various diffractive eigenstates i and k, for given separations b = bi — b2 
between the incoming protons. These are given by 

b) = ^ / n^,iy, bi, h2)ni,iY - y, bi, h2)d%d%5^'\h^ - b2 - b) (37) 

where F = Ins, see Fig. 6. Note that there is no integral^ over y. The convolution may be 
calculated at any rapidity y, leading to the same result. Given this effective Hk eikonal', we 
can calculate the cross sections (analogously to (5)- (7)). We obtain 



crtot(>^,b) = 2^ laiplflfcl 



(38) 



dcJoiiY, b) ^ _1_ 

dt 4:71 



i,k 



where t = —q^, and 



(39) 



aei(y,b) = J d% (i 

i,k 



,F,fe(y,b)/2 



(40) 



''Note that all the opacities in the absorptive exponents are multiplied by A = 0.25. Thus, the value of the 
product XvRR = 1 is not large. 

^Therefore it will be interesting and important to measure the single diffractive cross section, not only in the 
usual form dcgD/dM^, but also in the form dagD/drj, where rj denotes the position of the edge of the rapidity 
gap. This may be possible using the forward shower counters proposed in [28]. 

^The integral over y gives the multiplicity. 
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Figure 7: The symbolic diagram for (41) for low-mass dissociation of the 'beam' diffractive 
eigenstate i. 



4.3 Low-mass diffractive dissociation 



For low-mass excitation of the beam proton we obtain 

daei+sD{Y,h) 1 



dt 



Air 



El 1 2 1 1 2 1 I 



X 



X 



j d^'e-'^-"^' (l - e^'^'(^'b')/2)l U d^e'^-"" (l - e^^'^-^^'b)/^ 



(41) 



which has the symbolic structure shown in Fig. 7. 

Strictly speaking, we may need a different diagonalisation matrix a of (13) for the different 
t-channel exchanges. However, if the main difference between the diffractive eigenstates is due 
to the size and the impact parameter structure of the state, which is frozen for a fast hadron 
during the interaction, then it is justified to use the same eigenstates for any t-channel exchange, 



4.4 High-mass diffractive dissociation 

The expression for the high-mass excitation is more complicated. The cross section for beam 
particle diffractive dissociation (with the gap up to y) can be written using (5)-(7). Diffractive 
dissociation may be considered as the elastic scattering of intermediate parton c caused by its 
absorption on the target, which is described by the factor exp(— Ai7fc/2). 

Thus, in each impact parameter point b the cross section for single dissociation is propor- 
tional to (i) the elastic c—k cross section (1 — exp(— Af2fc(?/, b)/2))^; (ii) to the probability to find 
the intermediate parton c in the interval dy, that is Aexp(— Af2j/2 — Ai7fc/2); (iii) to the ampli- 
tude Vti of the parton c-beam interaction; (iv) to the gap survival factor S'^(b) = exp(— fi(y, b)) 
(y = Ins). The resulting cross section reads 

^^SD ^ /"q _ g-ACfe(y,bl,b2)/2^2^g-A^,{y-^/,bl,b2)/2-A^fc(J/,bl,b2)/2 ^ 

dy J 

X n,{Y - y, bi, b2)S'^(|bi - h,\)d%,d%2 , (42) 
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^g-A(Q,(y-y)+afe(j/))/2 

dashed lines 



Figure 8: The symbolic diagram for (42) for high-mass dissociation of the 'beam' diffractive 
eigenstate i. 



where bi (b2) are the coordinates in the impact parameter plane with respect to the beam 
(target) hadron. The normalisation factor is specified in (44). The gap survival probability^^ 

5f,(b)=exp(-F,,(b)). (43) 

The symbolic structure of (42), for high-mass single dissociation, is shown in Fig. 8. 

Accounting for the different Good- Walker eigenstates and the different states in t-channel 
we obtain 



2 ' - 2 



E 



a. 



|2 



E E T^tkiv. bi, ^2)Mk{v. bi, b2))5.,(|bi - b^l 



X 



X n^(F-|/,bi,b2)rf%rf262//?o, (44) 

where the parton density p"^ is defined by eq.(33) and, neglecting the secondary reggeon con- 
tribution, the elastic c — k amplitude 

Tt^{y, bi, b2) = (l - e-^^^('^'*'-''^)/2) . (45) 



For the secondary Reggeon, the real part may be large: note that tt/2((^)) ~ 0.2. To allow 
for this, and since, that besides the /2-trajectory with vacuum quantum numbers, there exists 
also tu, p and 02 exchange, we enlarge the contribution due to fi^ by increasing the values of 
the effective triple-Reggeon couplings gppR, guRp, Qhrr as compared to those coming from 
the absorptive opacities (31). So we use (45) for a = P2 or P3. On the other hand, for a = Pi 
or R we use, respectively. 



T^'iy, bi, b2) = V (1 - e-^"fe'(2''bi,b2)/2^ + r^^p (^i _ ^-xnR{yMM)vp^/2^j _ ^45) 



^"Strictly speaking, when calculating the gap survival probability in each particular case, we only have to 
account for the possibility of rescattering which produces secondaries within the gap interval. That is, in (43) 
we should not put the whole irreducible amplitude Fik{h), but, instead, part of it; since the contribution from 
the processes with a gap in the same (or a larger) rapidity interval does not change qualitatively the structure of 
the diffractive dissociation event. In the present computations we neglect this effect. This means that actually 
the gap survival probabilities, and the true cross sections of diffractive dissociation, should be a bit larger. 
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and 

7;f (y, bi, h,) = ^(l - e-^"'.(^.b„b.)/2)2 ^ ^^^^ _ ^_,nH(,,b„b.K,/2)2 _ (47) 

where = tpprD,^^ +fl^VRR. We take r^j^jp = rppR = 3, trrr = 9 to reproduce the available 
data in the CERN-ISR to Tevatron energy range. 

The slope of the diffractive dissociation cross section, Bsb = d[lia{das-D/ dM'^]/ dt at t = 0, 
can be calculated as the mean value of b| - the separation of the intermediate parton c from 
the target hadron 



SD 



k a 



X 



dM^ 



(48) 



4.5 Central exclusive production 



Central Exclusive Diffractive (CED) production of a system with mass squared = ^1^2^ with 
the large rapidity gaps either side, which is sometimes called the Double-Pomeron-Exchange 
(DPE) process, has a cross section given by 



d^id^2 



E 



i.k 



dXd%d%/(3l 



where 



Ef = \a,\'T^^{Y - y„ bi, h',)^ p^HY - y„ bi, b^ , 
E^k = \a,\'T^,{y2,h[M)^JpUy2MM) 



(49) 

(50) 
(51) 



are the probability amplitudes for elastic scattering of the intermediate parton c (c') on the 
beam (target) eigenstate i (k). The coordinates of parton c (c') are bi and ba {h[ and b2) 
with respect to the beam and the target proton respectively; that is, b'^ = bi — b2 + b2. The 
momentum fractions {^i — 1 — x^^i) of the incoming protons, transferred across the gaps, are 
^1 = e-^^-^'i) and ^2 = e-^^. The gap survival factor Sik{h = |bi - h'^\) is given by (43). 

The amplitude of the interaction of partons c and c', Q""'(yi, ^2, b2, b2)), is obtained by the 
solution of the evolution (29), which starts from the initial condition il"(y = ^2) = (^^^^ (bs — b2). 
That is, it starts from one parton at rapidity y2 placed at coordinate b2 in t-channel state a, 
and finishes at the point yi, h'2 in the state a'; note yi > y2- After the usual solution of (29,30), 
the evolution (29) was performed in the known "background" fields f)^, il" to account for the 
absorption of intermediate partons. 
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da,,/dt (mb/GeV^) 




-t (GeV^) 

Figure 9: The t dependence of the elastic pp cross section. The dashed and dotted hnes are the 
contributions from the elastic scattering of the largest size (i = 1) and the smallest size {i = 3) 
components. 

5 Description of the data and predictions for the LHC 

Clearly, the number of parameters in our model is too large to perform a straightforward fit of 
the data. Instead, we fix the majority of the parameters at reasonable values and demonstrate 
that such a model can reproduce all the features of the available data on diffractive cross 
sections, (Jtot, daei/dt, crgg"^, dasi)/dM^. 

Before we give the values of the parameters, we show in Fig. 9 the quality of the description 
of the data for the elastic differential cross section. We also present in Fig. 9 the prediction for 
differential elastic cross section at the LHC energy ^/s — 14 TeV. Recall that we are using a 
three-channel eikonal. That is i,k — 1,2,3. It is interesting to note that the contribution to 
the cross section arising from the scattering of the two large-size eigenstates, (i = 1) x (i = 1), 
already has a diffractive dip at —t = 0.2 GcV^. However, after the contributions from all 
possible combinations i x k are summed up, the prediction has no dip up to —t = 0.5 GeV^. 
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Note, also, that the Pomeron and secondary Reggeon couphngs to the proton were taken to 
have the forms 

(3(t) = Poe'^'/il - t/d^f, m = ^Re^-K (52) 

The values of the parameters that we use are 

d2 = 0.15 GeV-^ di = 1.5 GeV^ dn = 1 GeV'^. (53) 

The non-zero value of ^2 is simply to provide good convergence and accuracy of the Fourier 
transform. The parameter di controls the t behaviour of the elastic cross section, while du is 
responsible for the t slope of diffractive dissociation at relatively low y = —In^, where the cross 
section is dominated by the RRP triple-Reggeon term. The relative size of this contribution, 
as compared to that due to PPP, was tuned by choosing vpr = 1.8 and r/jRp = 3. In order 
to describe the data, the couplings were found to be /Sq = 33 mb and jS'^ = 8 mb. Since wc 
choose a relatively simple t dependence for the Reggeon-proton couplings (3{t), our model is 
only reliable over a restricted t interval, —t < 0.5 GeV^. Note that in this domain, the real 
part of Pomeron exchange, and a possible Odderon exchange contribution, would give only very 
small effects. 

To describe the high energy behavior of the total cross section, we take A" = 0.3 for each of 
the three components of the Pomeron. These Pomeron intercepts are consistent with resummed 
NLL BFKL, which gives coq ~ 0.3 practically independent of the scale kt [29]. The slopes of the 
Pomeron trajectories are driven by the transverse momentum associated with the particular 
component a. In fact, we have a' oc l/k^. We find the data require a'p_^ = 0.05 GeV~^ for 
the large-size Pomeron component, so we put a'p^ = 0.05/9 GeV^^ for the second component 
and a'p^ = for the smallest-size component. For the secondary Reggeon trajectory we take 
a'ji = 0.9 GeV~^, and ^^(O) = 0.6. The 'bare' value is a bit larger than |, since the final 
effective intercept is reduceed by the absorptive corrections included in the evolution equation. 
The description of the total cross section data are shown in Fig. 10(a). The screening corrections 
arising from the 'enhanced' multi-Pomeron diagrams, that is from the high-mass dissociation, 
slow down the growth of the cross section with energy. Thus, the model predicts a relatively 
low total cross section at the LHC - (Ttot(LHC) ~ 90 mb^^. 

5.1 Low-mass dissociation 

Recall that the couplings of the Good- Walker eigenstates i were specified by Pi(t) = ji(3{'yit). 
The values 71 = 1.80, 72 = 0.82 and 73 = 0.38 were chosen so as to reproduce the low-mass 
dissociation cross section crgg"^ = 2 mb at the CERN-ISR energy^^ [27]. 

-^-"^This value is also predicted by other models of 'soft' interactions which include absorptive effects [30, 31]. 
Although, here, we use a three-channel eikonal model, practically the same results, and the same quality of 
the description, is obtained using a two-channel eikonal, that is only two eigenstates (see also the discussion 
in [22]). 
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5.2 High- mass dissociation 

The value of the parameter A, which controls the cross section of high-mass dissociation in the 
small ^ (that is, large y) region, was found to be A = 0.25. The dependence of the cross section 
for high-mass dissociation, ^(Pa/dtd^, on ^ = M^/s is compared with the Tevatron CDF data 
[32, 33] in Fig. 11. Recall that in our model we have not included pion exchange. Since the 
ttttP term is essential at large ^, we have included the corresponding contribution using the 
parameters obtained in [7]. The results without the ttttP term arc shown by the dashed lines. 
We also show in Fig. 11(a) (by the dotted line at small ^) the prediction for the LHC energy. 

Above, we introduced the different types of data mentioning the parameters that they 
mainly constrain. Of course, in practice, these parameter values are used to describe all the 
'soft' data simultaneously. 

The energy behaviour of the cross sections are shown in Table 1 and Fig. 10. Fig. 10 also 
shows the energy behaviour of the multiplicities of the secondaries produced by the i-channel 
Pomeron components of different sizes; we will discuss the multiplicity distributions in some 
detail in Section 6. 
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tot 
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73.7 


16.4 


4.1 


9.7 


13.8 


14 


91.7 


21.5 


4.9 


14.1 


19.0 


100 


108.0 


26.2 


5.6 


18.6 


24.2 



Table 1: Cross sections (in mb) versus collider energy (in TeV). 

The values of quoted in Table 1 look, at first sight, too large, when compared with 
the value 9.46 ± 0.44 mb given by CDF [32]. However the CDF value does not include the 
secondary Reggeon (RRP) contribution, denoted as a 'non-diffractive' component of 2.6 ± 0.4 
mb. Moreover, the trigger used to select the diffractive dissociation events rejects part of the 
low-mass proton excitations. Taking these absences into account, there is no contradiction 
between the model prediction and the CDF data. Furthermore, note that in the region where 
the CDF detector efficiency and resolution are good, our model gives an excellent description 
of the measured data, sec Fig. 11. 

It is interesting to note, that after tuning the parameters to describe all the available 'soft' 
data, the model satisfies the Finite Energy Sum Rules [34] to good accuracy^"^. Indeed, we can 
switch off the low-mass dissociation, putting the same couplings for each diffractive eigcnstate 
(71 = 72 = 73 = 1), and replace the low-mass excitations for Ay < 1.5 by the triple- and multi- 
Regge contributions-*^^. Keeping all the other parameters as before, we then obtain (Jtot =73 mb 

"^^We thank Alan White for discussions. 

^^Recall that in the basic model we introduced a threshold Ay = 1.5; we started the evolution (29), (30) at 
y = 1.5 in order not to generate low-mass dissociation via triple- and multi-Regee contributions and to avoid 
double counting. 
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(93 mb) and f7g^=13.6 mb (20.1 mb) for the Tevatron (LHC) energies, These values are close 
to those in Table 1. 

In principle, it is straightforward, although computer intensive, to use the model to calculate 
the cross section for double dissociation, ctdd- We do not show the values here. The values will 
be similar to those in Table 2 of Ref. [22]; footnote 25 of that paper shows that the values of 
(Tdd are in excellent agreement with the Tevatron data. 

5.3 Central exclusive production 

The cross sections for the Central Exclusive Diffractivc (CED) production at the LHC energy 
■\/s = 14 TeV are shown in Fig. 12 for those ^ intervals which can be studied by the TOTEM 
and FP420 detectors. Here we mean the soft CED production of a state with the mass given by 
= ^i^2S separated from the incoming protons by two large rapidity gaps. The calculation is 
described in Section 4.5. The cross section integrated over the 0.002 < < 0.2 (0.02 < < 0.2) 
intervals (for both and ^2) is predicted^^ to be 53 (16) fih. The major contribution comes 
from pure soft interactions. For ^ < 0.02, more than 80% of the cross section is due to the 
large size component of the Pomeron; and more than half for larger values of ^. Note that, in 
the CED calculations, we did not include the vrvrP contribution. Thus, actually, the expected 
cross section will be larger for ^ values that are not too small, see Fig. 11. 

Note that the resulting CED cross section is about twice larger than that expected from 
the naive factorization formula 

d^idC.2 o-tot d^i d^2 

This is due to the fact that each single dissociation contains a gap survival factor S%{h) and 
therefore the r.h.s. of (54) is proportional to {S'^Y while the l.h.s. contains S"^ once only. On the 
other hand, for Central Exclusive Production the typical values of the impact parameter h are 
smaller; so we have a smaller gap survival factor (5'^(6)) . This partly compensates the additional 
power of 5"^ in the r.h.s. of (54), and as a result the violation of the factorization shown in (54) 
is not so strong. We will discuss rapidity gap survival in central exclusive production in detail 
in the following paper [15]. 



6 Multiparticle inclusive production 

As we have a detailed model for high energy soft processes, it would appear to be possible to 
predict the multiplicity distribution at the LHC. However, although some general features can 

^^To speed up the computation we neglect the small non-zero value of a' in the calculation of the amplitude 
Qaa Then there is no diffusion in impact parameter space and the integral over h'2 disapears; since h'2 = b2. 

However, to correct the final result w{^ smear o\ii the resiilting amplitude, which allows for the larger gap survival 
probability at larger h. In this way we retain reasonable (~ 20%) accuracy of the computations. 
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be predicted, it is not so simple to make a quantitative prediction. We address the problem 
below. 

Recall that in the evolution equations for the amplitude, given in (29), (30), we include the 
absorptive factor exp(— Q/2), and not exp(— Q). That is we work with the forward amplitude 
ImT(6) — 1 — e~^/^, which at each step of the evolution (in rapidity y) includes all possible 
processes - both elastic and inelastic interactions with cross sections cTci{b) = (1 — e~^/^)^ and 
cr-meiib) = 1 - e~^; where atot{b) = 21mT{b) = aei{b) + ainei{b), see eqs. (5)- (7). 

As usual, inelastic processes include both single-ladder exchange, as well as multiple inter- 
actions with a larger density of secondary partons. The situation is similar to the rescattering 
of a fast hadron in a heavy nucleus. That is, in such an eikonal approach the probability, WN{b), 
of events with parton multiplicity N times larger than that in a single ladder, is given by 

= ^e-^. (55) 

In the multi-channel case the opacity Q should be replaced by Fn-. Unfortunately, we cannot use 
this probability Wn literally to describe the multiplicity distributions of secondary hadrons. In 
particular a non-negligible fraction of the final hadrons may be produced via the fragmentation 
of minijets. These processes are beyond the 'pure soft' approach used in the present paper. 
Therefore, below, we discuss the multiplicity distribution only at the partonic level. 

The mean number of the (t-channel) ladders of the type a produced in the collision of i and 
k Good- Walker eigenstates can be calculated as 

NtM = o2 / bi, b2)Jlf (y - y, bi, b2)d'6id^625('nbi - b2 - b) . (56) 

where (7ifc(b) = 2[1 — exp(— Fjfc(b)/2)] and a — Pi, P2, P3, R. Recall that Pi, P2 and P3 are the 
large, intermediate and small size components of the Pomeron respectively. After averaging 
over the impact parameter b and the diffractive eigenstates i, k of the incoming protons, we 
obtain 

Kube = -^J2\ai\'\a,\' [ni{y,hi,h2)nt{Y-y,hi,h2)d%id%2. (57) 

This quantity may be considered as the mean number of colour tubes of type a produced in the 
proton-proton interaction. Note that the value of A^tube does not depend^^ on the rapidity y. 

To obtain the number of partons created by the ladder 'a' at rapidity y, we have to include 
the parton density p"(y) of (33) in the numerator of (57). That is 

^p^rton = E I niiV: bi, h2)pt,n^{Y - y, hi, h2)d%d%2 ■ (58) 



The results are shown in Fig. 10(d). The main growth in multiplicity, as we go from Tevatron to 
LHC energies, is due to the small size ('QCD') Pomeron component, which produces particles 
with typically pt ^ 5 GeV. There is essentially no growth in multiplicity at small pt- This 
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simply confirms the trend that has been observed through the CERN-ISR to Tevatron energy 
range, see the data points in Fig. 13. 

In other words, starting with the same intercepts (A = 0.3) the large size component contri- 
bution after the absorptive correction becomes practically fiat, while the small size contribution, 
which is much less affected by the absorption, continues to grow with energy. As mentioned 
above, such a behaviour is consistent with the experiment (see Fig. 13) where the density of a 
low kt secondaries is practically saturated while probability to produce a hadron with a large 
(say, more then 5 GeV) transverse momentum grows with the initial energy. 

Unfortunately we cannot identify each 'parton' with a pion. At a large kf it is probable 
that, after the hadronisation, a parton forms a gluon jet of pions. However, at low kt this is 
not evident; it is hard to say what a gluon 'jet' becomes at a low kt. Nevertheless, in order 
to compare the model with the data we assume, that after hadronisation, each parton from 
Pomeron components Pi, P2 and P3 gives three charged pions with p^j ~ 0.5, 1.5 and 5 GeV 
respectively. In this way we may estimate the inclusive cross section, at the three values of pt, 
using 

Eda'' _ 3crtotA^partoii /rqv 
d^p npl 

where Utot allows for normalisation and 1/pf accounts for the size of the phase space occupied 
by the particles from component i. 

To obtain a qualitative feel for the expected behaviour, we show our predictions at the 
Tevatron and LHC energies in Fig. 13, where the horizontal lines indicate the typical pt interval 
associated with each Pomeron component. Recall, that in our 'soft' model, we never use 
the value of the Pomeron kt explicitly. The characteristic parameters actually used in the 
computations are the ratios kl ^/k^ ^j^^. The horizontal lines reflect the pt intervals covered by 
the various components arising from the scale choices. Some features of Fig. 13 are clear. First, 
although there is some freedom in assigning the overall scale, nevertheless, it appears that the 
scale choice made in the figure agrees satisfactorily with the Tevatron data [35]. Second, as 
compared to the Tevatron, the LHC distribution is more enhanced at large pt. The enhancement 
is a factor of 2.6 for the 'QCD' small-size component of the Pomeron, whereas it is only 1.25 
for the 'soft' Pomeron component. 



7 Summary 

New triple- Regge analyses [7, 36], which include absorptive effects, found that the triple- 
Pomeron coupling is rather large {g^p — Xqn with A > 0.2). Thus, in order to obtain reli- 
able predictions for diffractive processes at the LHC, it is necessary to have a model of 'soft' 
high-energy processes which includes multi-Pomeron interactions. 

^^This was checked by straightforward computation. 
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Here we have presented such a model, tuned to the existing 'soft' data, which, in principle, 
is capable of predicting the basic features of high-energy soft pp interactions. The absorption 
of intermediate partons is described by conventional exp(— AO)-type factors. This corresponds 
to a coupling — nm\^^'^~'^gN of the n ^ m Pomeron vertices. 

Briefly, the model has multi-components in both the s- and t-channels. The former are based 
on a three-channel eikonal approach, together with the inclusion of multi-Pomeron diagrams, 
so that both low- and high-mass diffractivc dissociation are well described. Predictions for the 
LHC are given. A novel feature of the model is the inclusion of different t-channcl exchanges, 
which allows for small-, intermediate- and large-size components of the exchanged Pomeron, 
each with a bare intercept A = ap{0) — 1 — 0.3. For the large-size component, the slope of 
the trajectory is a'p — 0.05 GeV~^. The large-size Pomeron component is heavily screened by 
the effect of 'enhanced' multi-Pomeron diagrams, associated with high-mass dissociation. This 
leads, among other things, to the effective "saturation" of the low pt particle density, and to a 
slow growth of the total cross section. Indeed, the model predicts a relatively low total cross 
section at the LHC - o'tot(LHC) ^ 90 mb. On the other hand, the small-size component of 
the Pomeron is weakly screened, leading to an anticipated growth of the particle multiplicity 
at large pt (~ 5 GeV) at the LHC. Thus the model has the possibility to embody a smooth 
matching of the perturbative QCD Pomeron to the 'soft' Pomeron. 

We emphasized that a reliable model of soft interactions is essential in order to predict the 
rates of diffractive processes at the LHC. In particular, we used the model to calculate the 
rapidity gap survival factors, including the effects of both eikonal and enhanced rescattering. 
This is the subject of the following paper [15]. 
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Figure 10: The energy dependence of the total (a), elastic and diffractive dissociation (b) pp 
cross sections and the cross sections of dissociation to a fixed = state (c); (d) the parton 
multiphcity (solid lines) and the number of 'colour tubes' (dashed) produced by the Pomeron 
components of different size. 
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Figure 11: The model description of the data for the cross section for high- mass dissociation 
versus ^ for -t = 0.05 GeV^ at = 1800 GeV and 546 GeV [32, 33]. The dashed hues are 
the predictions without the vrvrP contribution. The dotted curve at small ^ is the prediction 
for the LHC. 
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Figure 12: Sample predictions for Central Exclusive Diffractive production at the LHC. The 
^j's are the momentum fractions of the incoming protons transferred across the rapidity gaps 
on either side of the centrally produced system of mass M — y/^i^. 



27 



> 

O 



o 



E 

o 



CL 

to 

■D 



b 
■a 

UJ 



10 



-24 



10 



-26 = 



10 



-28 



10 



-30 : 




10 



-32 



10 



-34 



T 1— 1 1 — 

O CDF 1 800 GeV 
AUA1 546 GeV 
□ BS 53 GeV 



27 GeV ^ 



10 



-36 







LHG 



□ ^ — 



^Tevatron 



0.5*(h'+h') 



4 6 
Pt (GeV/c) 



8 



10 



Figure 13: The plot is from Ref. [35]. The horizontal lines, which are superimposed, are our 
model predictions at the Tevatron and LHC energies; the three pt ranges correspond to the 
large-, intermediate- and small-size components of the Pomeron. 



28 



